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Abstract. When can a model of a physical system be regarded as com- 
putable? We provide the definition of a computable physical model to answer 
this question. The connection between our definition and Kreisel's notion of a 
mechanistic theory is discussed, and several examples of computable physical 
models are given, including models which feature discrete motion, a model 
which features non-discrete continuous motion, and probabilistic models such 
as radioactive decay. We show how computable physical models on effective 
topological spaces can be formulated using the theory of type-two effectiv- 
ity (TTE). Various common operations on computable physical models are 
described, such as the operation of coarse-graining and the formation of sta- 
tistical ensembles. The definition of a computable physical model also allows 
for a precise formalization of the computable universe hypothesis — the claim 
that all the laws of physics are computable. 



1. Introduction 

A common way to formalize the concept of a physical model is to identify the 
states of the system being modeled with the members of some set S, and to identify 
each observable quantity of the system with a function from S to the real numbersu 
For example, a simple model of planetary motion, with the Earth moving in a 
circular orbit and traveling at a uniform speed, is the following. 

Model 1.1 (Simple Planetary Motion). Let S be the set of all pairs of real numbers 
{t,a) such that a = 360 (t — [ij), where \ t\ denotes the largest integer less than or 
equal to t. The angular position of the Earth, measured in degrees, is given by 
the function a{t,a) = a. The time, measured in years, is given by the function 
T{t, a) — t. 

If we wish, for example, to compute the position of the Earth after 2.25 years, 
we ask: "For which states (t,a) does T{t,a) = 2.25?" There is only one such 
state, namely (2.25,90). Therefore, the position of the Earth after 2.25 years is 
a(2.25, 90) — 90 degrees. We say that the model is faithful if and only if the values 
of the observable quantities in the model match the values that are physically 
observed. 

Church and Turing hypothesized that the functions which are effectively com- 
putable by humans are exactly the recursive functionsll There have been several 
attempts [HI [29l 13 [28l [21] to extend the Church- Turing thesis to physics, hy- 
pothesizing that the laws of physics are, in some sense, computable. But given an 
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arbitrary physical model, it has not been clear exactly how one determines whether 
or not that model is to be regarded as computable. To date, the best attempt at 
providing such a definition has been Kreisel's notion of a mechanistic theory [llj . 
Kreisel suggested the following. 

Kreisel's Criterion. The predictions of a physical model are to be regarded as 
computable if and only if every real number which is observable according to the 
model is recursive relative to the data uniformly. 

But many seemingly innocuous models have failed to satisfy Kreisel's criterion. For 
example, the simple model of planetary motion (Model II. 1|) fails because given a 
real number representing the time t, there is no effectively computable procedure 
which determines the corresponding angle a when a is near the discontinuity at 360 
degrees. Models which intuitively seem to have computable predictions often fail 
to satisfy Kreisel's criterion because discontinuities in their formalisms prevent the 
models' predictions from being regarded as computable, despite the fact that there 
are no discontinuities in the actual physical phenomena being modeled [TB] . 

Rather than using Kreisel's criterion to prove that the predictions of established 
models are computable, an alternate approach is to supply a restrictive formalism 
which guarantees that the predictions of models expressible in that formalism are 
computable. This has been the approach taken in references [HI [29l \8\ [28] . But 
difficulties have been encountered expressing important established models in these 
formalisms. For example, Rosen [18] was unable to describe radioactive decay in 
the formalism that he had proposed, and work is ongoing to describe established 
physical models in other computable formalisms. 

It is the goal of this paper to provide a general formalism for describing physical 
models whose predictions are computable, and to show that the computable for- 
malisms studied by previous authors are special cases of our general formalism. In 
particular, we show in Section [TS] that among the members of a large class of phys- 
ical models, each physical model satisfying Kreisel's criterion has a corresponding 
model in our formalism. We also avoid some of the difficulties which, for example, 
prevented the simple model of planetary motion (Model II. ip from being regarded 
as computable, as will be seen in Section S] Our approach also avoids the difficulty 
that Rosen encountered with radioactive decay, as will be seen in Section [B] 

2. Computable Physical Models 

The central problem is that physical models use real numbers to represent the 
values of observable quantities, but that recursive functions are functions of non- 
negative integers, not functions of real numbers. To show that a model is com- 
putable, the model must somehow be expressed using recursive functions. Careful 
consideration of this problem, however, reveals that the real numbers are not actu- 
ally necessary in physical models. Non-negative integers suffice for the representa- 
tion of observable quantities because numbers measured in laboratory experiments 
necessarily have only finitely many digits of precision. For example, measurements 
of distances with a measuring stick will always be non-negative integer multiples of 
the smallest division on the measuring stick. So, we suffer no loss of generality by 
restricting the values of all observable quantities to be expressed as non-negative 
integers — the restriction only forces us to make the methods of error analysis, which 
were tacitly assumed when dealing with real numbers, an explicit part of each 
model. 
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Non-negative integers are not only sufficient for tire description of direct physical 
measurements, but are also sufficient for encoding more complex data structures — 
allowing us to define recursive functions on those data structures. For example, a 
pair of two non-negative integers x and y can be encoded as a single non-negative 
integer (x, y) using Cantor's pairing function 

{x, y) = ^(x^ + 2xy + y'^ + ix + y) 

A pair (a;, y) of non- negative integers will also be called a length two sequence of non- 
negative integers. A triple (or equivalently, length three sequence) of non-negative 
integers x, y, and z can be encoded as {{x, y), z), and so on. We write {x, y, z) as an 
abbreviation for {{x,y),z). An integer i can be encoded as a non- negative integer 
C(i) using the formula 

f-2f-l ifi<0 

C(«) = S„. ..... 

\2i if I > 

And a rational number | in lowest-terms with 6 > can be encoded as a non- 
negative integer /3(f) using the formula 

p(^) ^ ^((sgna)2C(''i-''i)3«("^-''^)5^(''='-''^)7^("*-'"'hl';(''=-''^) • • • ) 

where a = (sgna)2°i3''^5°^7°*'ll''^ • • • is the prime factorization of the integer a, 
and similarly for h. We write [q ; r) as an abbreviation for the pair of rational 
numbers (^p{q) , p{r)'^ . 

Historically, authors who have wished to restrict themselves to physical mod- 
els whose predictions are computable have chosen from among a handful of for- 
malisms. For example, Zuse [52] and Fredkin |B] have formalized their models as 
cellular automata, with each cell of an automaton representing a discrete unit of 
space and each step of computation in the automaton representing a discrete unit 
of time. Wolfram [35] has formalized his models in a variety of computational sys- 
tems, including cellular automata, but has favored network systems for a model 
of fundamental physics. In each of these cases, the states of a physical system 
are represented by the states of a computational system (for example, a cellular 
automaton or a network system) which can be encoded as non-negative integers 
using the techniques just described. The resulting set of non-negative integers is a 
recursive set, and the observable quantities of the system are recursive functions of 
the members of that set. This immediately suggests the following definition. 

Definition 2.1. A computable physical model of a system is a recursive set S of 
states with a total recursive function (j) for each observable quantity of the system. 
(/)(s) is the value of that observable quantity when the system is in state s. 

So, in a computable physical model the set 5' is a set of non- negative integers, and 
each observable quantity is a function from non-negative integers to non-negative 
integers. The models considered by Zuse, Fredkin, and Wolfram are necessarily 
special sorts of computable physical models, and the set of all computable physical 
models is a proper subset of all physical models. In order to avoid all ambiguity, we 
insist that observable quantities be defined operationally [2] in computable physical 
models, so that, for example, if there were an observable quantity corresponding to 
time, then that observable quantity would be the time as measured with a specific 
conventionally-chosen clock in a specific conventionally-chosen reference frame. 
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An immediate consequence of the definition of a computable pliysical model is 
that we can give a precise formal counterpart to the informal claim that all the 
laws of physics are computable. 

Computable Universe Hypothesis. The universe has a recursive set of states 
U . For each observable quantity, there is a total recursive function <f>. (j)(s) is the 
value of that observable quantity when the universe is in state s. 

By a distinguishable system, we mean any system for which there is an observable 
quantity </> such that 0(s) = 1 when the system exists in the universe, and such that 
(j){s) = otherwise. For example, if the system being studied is the orbit of the 
Earth, then (f){s) = when state s corresponds to a time before the formation of the 
Earth, and (j){s) — 1 when the Earth exists and is orbiting the Sun. Note that the 
set of states s inU for which 0(s) = 1 is itself a recursive set whenever U and (p are 
recursive. So, the computable universe hypothesis implies that computable physical 
models are sufficient for modeling any distinguishable system in the universe — the 
set of states of that distinguishable system is the set of all members s of U for 
which (/)(s) = 1, and the observable quantities of the distinguishable system are 
necessarily a subset of the observable quantities of the universe. 

3. Discrete Planetary Motion 

As a first example of a computable physical model, consider the following model 
of planetary motion. 

Model 3.1 (Discrete Planetary Motion). Let S be the set of all pairs {{r ; s) , {p ; q)') 
such that 

* ^ p = 360(r-[rj) 



10 100 

i + 1 1 



q = 360(3- [s\) 



10 100 

for some integer i between —20000 and 20000. The angular position of the Earth, 
represented as a range of angles measured in degrees, is given by the function 
a((r ; s), (p ; (7)) = {p',q)- The time interval, measured in years, is given by the 
function T((r ; s), (p ; g)) = (r ; s). 

This is a discrete model. That is, the position of the Earth in its orbit is not an ex- 
act real number, such as 90 degrees, but is instead an interval, such as (68.4 ; 111.6) 
representing a range of angles between 68.4 degrees and 111.6 degreeso Similarly, 
time is measured in discrete intervals of length 0.12 years. The earliest time in- 
terval in the model is near the year —2000 and the latest time interval is near 
the year 2000. Moreover, this model is faithful — it is in exact agreement with all 
observations. 

There are ten possible measurements for the angular position of the Earth in 
Model O 



(32.4; 75.6) (68.4; 111.6) (104.4 ; 147.6) (140.4 ; 183.6) (176.4 ; 219.6) 
(212.4; 255.6) (248.4 ; 291.6) (284.4 ; 327.6) (320.4; 3.6) (356.4; 39.6) 



■^We use decimal numbers to represent exact rational numbers. For example, 68.4 is to be 
understood as an abbreviation for 



THE COMPUTABLE UNIVERSE HYPOTHESIS 



5 



These are the intervals obtained by dividing the 360 degrees of the circle into 
ten equal intervals of 36 degrees each, then extending each interval by exactly 3.6 
degrees on both sides, bringing the total length of each interval to 43.2 degrees. 
Therefore, consecutive intervals overlap by 7.2 degrees (there is also overlap in 
consecutive time intervals), and this serves an important purpose. The Earth's 
orbit is not, in reality, a perfect circle, and the Earth does not spend an equal 
amount of time in each of the intervals. But because the eccentricity of the Earth's 
orbit contributes to, at most, only about a 2 degree deviation [7 from the simple 
model of planetary motion (Model [TTT|) . the overlap of these intervals is more than 
adequate to conceal evidence of the eccentricity, ensuring that this discrete model 
is faithful. Also note that the overlap is a realistic feature of all known instruments 
which measure angles, since each such instrument has only a limited accuracy. If 
angles are measured with a protractor, for example, the accuracy might be limited 
by the thickness of the lines painted on the protractor, which divide one reading 
from another. For example, if the lines are 7.2 degrees thick, then it might not be 
possible to distinguish a reading of (32.4; 75.6) from a reading of (68.4; 111.6) if 
the quantity being measured is somewhere on that line (that is, if the quantity is 
somewhere between 68.4 and 75.6 degrees). The accuracy of measuring instruments 
is discussed in greater detail in Section [SI 



4. Non-Discrete Continuous Planetary Motion 

Many commonly-studied computable physical models are discrete, but non- 
discrete continuous models are also possible. For example, a non-discrete continuous 
computable physical model of planetary motion is the following. 

Model 4.1 (Non-Discrete Continuous Planetary Motion). Let S be the set of all 

pairs ((r ; s), (p ; g)) such that 

* p = 360(r- [rj) 



10" 10"+i 
i + l 1 



g = 360(s- [sj) 



10" 10"+i 

for some integer i and some positive integer n. The angular position of the Earth, 
represented as a range of angles measured in degrees, is given by the function 
a((r ; s), (p ; q)) = (p;<z)- The time interval, measured in years, is given by the 
function T((r ; s), (p ; g)) = (r ; s). 

Like the discrete model, angular position and time are measured in intervals, 
but in this case the intervals are not all the same length. In particular, there 
are arbitrarily small intervals for the observable quantities of position and time, 
meaning that these quantities may be measured to arbitrary precision. This feature 
of Model 14.11 allows us to speak about real-valued positions and times, despite 
the fact that the values of observable quantities in the model are all non-negative 
integers, not real numbers. 

This is because a real number is not the result of a single measurement, but is 
instead the limit of a potentially-infinite sequence of measurements. Suppose, for 
example, that we wish to measure the circumference of a circle whose diameter is 
exactly one meter. Measured with unmarked metersticks, we measure the circum- 
ference to be 3 meters. If the sticks are marked with millimeters, then we measure 



6 



MATTHEW P. SZUDZIK 



the circumference to be about 3.141 meters. And if they are marked with microm- 
eters, then we measure a circumference of about 3.141592 meters. If we continue 
this process indefinitely with increasingly precise measuring instruments, then in 
the infinite limit we approach the real number tt. 

More formally, for each real number x there is an infinite sequence of nested 
intervals (oq ; bg), (ai ; 6i), (a2 ; 62), . . . that converges to x. Given such a sequence, 
the function (p such that 4'{n) = (a„ ; 6„) for each non-negative integer n is said to 
be an oracle for x. Note that there is more than one distinct sequence of nested 
intervals converging to x, and therefore more than one oracle for each x. Of par- 
ticular importance is the standard decimal oracle Ox for the real number x. By 
definition, Ox{n) — (a„ where 

_ [lO"+ia;J c _ [lO"+ia;J + 1 c 



10"+i 10"+i 10"+i 10"+i 

for each non-negative integer n, and where the accuracy factor c is a positive 
rational number constant. We say that a; is a recursive real number if and only 
if Ox is a recursive function. Note that not all real numbers are recursive [22] . 

Now, returning to Model 14.11 suppose that we are asked to find the position 
of the Earth at some real-valued time t. Suppose further that we are given the 
oracle Ot with accuracy factor c= j^. Note that as we increase n, the values 0t{n) 
are increasingly precise measurements of the time t in Model 14.11 Therefore, for 
each n there is some state {ot{n),{pn;qn)) in the set S of Model |4?T] Because 
5' is a recursive set, and because there is exactly one state corresponding to each 
time measurement, the function e such that e(n) = (p„ ; (7„) is a recursive function 
relative to the oracle o*. In fact, if Ot{n) = (r„ ; s„), then 



e{n) = (360 (r„ - [r„J) ;360(s„ ~ 



for each non-negative integer n. And since the sequence of intervals (ro;so)j 
(ri;si), (r2;s2), ... converges to t, it immediately follows that the sequence of 
intervals e(0), e(l), e(2), ... converges to a = 360(t — [ij) whenever t is not an 
integer. In other words, e is an oracle for the angular position a. 
But in the case that t is an integer, 

^ ^ V 10"+i 10"+!/ 

for all non-negative integers n, and (356.4 ; 39.6), (359.64 ; 3.96), (359.964 ; 0.396), . . . 
is the resulting sequence. In the standard topology of the real numbers an interval 
{x ; y) should have x < y, so the question of whether or not this sequence converges 
to a point a in that standard topology cannot be meaningfully answered. But if 
we are willing to abandon the standard topology of the real numbers, then we may 
conventionally define this sequence to converge to a = 0. In fact, this definition 
is tantamount to establishing the topology of a circle of circumference 360 for all 
angles aQ Of course, this definition is justified since the readings after 360 on a 
measuring instrument for angles are identified with those readings after 0. In other 
words, angles really do lie in a circle. 



basis for this topology is represented by tiic set of all possible angle measurements. In 
particular, if x < y then (x ; y) represents the set of all real numbers a such that x < a < y, and 
if X > y then (x ; y) represents the set of all real numbers a such that 0<a<yorx<a< 360. 
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So, given the oracle ot for a real- valued time t, Model I4?l1 allows us to compute 
an oracle e for the angular position a of the Earth at that time. These predictions 
are in complete agreement with the predictions of the simple model of planetary 
motion fModel [O]) . In fact, imposing the appropriate topology on the space of 
angles a, the mapping from t to a in Model |4T1 is continuous. The same mapping is 
discontinuous in the standard topology of the real numbers, which leads Kreisel's 
criterion to fail for Model 11.11 The formulation of computable physical models on 
effective topological spaces is discussed in greater detail in Sections [TT] through [T5l 



5. Coarse-Graining 

Observable quantities in computable physical models are defined operationally. 
This means that each observable quantity is defined so as to correspond to a specific 
physical operation, such as the operation of comparing a length to the markings on 
a meterstick (where the meterstick itself is constructed according to a prescribed 
operation) . This is problematic for the non-discrete continuous model of planetary 
motion (Model 14. 1|) because, for example, arbitrary precision angle measurements 
are made with a single observable quantity in the model. That is, to assert that a 
model such as Model ITT] is faithful, one must assert that there exists an operation 
which is capable of measuring angles to arbitrary precision. It is not known whether 
or not such an operation actually exists. And although the point is somewhat moot, 
since Model 14.11 is clearly not faithful, it raises the question of whether this is an 
accidental feature of Model lTTTl or whether it is a feature common to all non-discrete 
continuous computable physical models. 

A more practical alternative to Model 14.11 might introduce an infinite sequence 
of observable quantities ai, a2, 03, • ■ • , each with finitely many digits of precision, 
and each more precise than its predecessor in the sequence. In this case, given 
a state s, the values ai{s), a2(s), a3(s), ... would form a sequence of intervals 
converging to a real number representing the angular position of the Earth in that 
state. But a computable physical model has only countably many states, and there 
are uncountably many real numbers in the interval (0 ; 360). Therefore, there must 
be some real number position in the interval (0 ; 360) that the Earth never attainsQ 
That is, a computable physical model of this alternative form is not continuous in 
the intended topology. 

Rather than considering arbitrary precision measurements, let us introduce just 
one additional level of precision into the discrete model of planetary motion (Model 
lO) . 



In particular, this is a real number constructed by diagonalizing over those real numbers which 
are associated with each of the countably many states. 
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Pi 


= 360 (ri - 


kij) 


qi 


360(si - 


bij) 


P2 


= 360 (ra - 


Vr2\) 


12 


= 360 (s2 - 


LS2J) 



Model 5.1. Let S be the set of all quadruples ((ri ; si), (r2 ; S2), {pi ', qi), {p2 92)) 
such that 

i 1 
''^ ~ 10 ~ 102 
+ l 1 

" ~W 102 

^ J L 

102 103 
102 103 

/or some integers i and j with lOi < j < lOi + 9. The angular position of the Earth, 
represented as a range of angles measured in degrees with a low-precision measuring 
instrument, is given by the function 

ai{{ri;si),{r2;s2),{pi;qi),{p2;q2)) = {pi;qi) 

The angular position of the Earth, represented as a range of angles measured in 
degrees with a high-precision measuring instrument, is given by the function 

(^2{{ri ; Si), (r2 ; S2), {pi ; qi), {p2 ; 92)) = (p2 ; 92) 

The time interval, measured in years by a low-precision measuring instrument, is 
given by the function 



((ri ; si), (r2 ; S2), [pi ; qi), {p2 ; 92)) = {ri ; si) 



The time interval, measured in years by a high-precision measuring instrument, is 
given by the function 

T2{{ri ; si), (r2 ; S2), {pi ; qi), {p2 ; '72)) = {r2 ; S2) 

Note that if the high-precision observable quantities a2 and T2 are ignored, then 
the predictions of Model 15.11 agree exactly with the predictions of Model 13. Ill The 
process of removing observable quantities from a model to obtain a new model 
with fewer observable quantities is called coarse- graining. But while Model [3TT] is 
faithful, Model 15.11 is not faithful — physical measurements do not agree with the 
values of the observable quantities a2 and T2 because the orbit of the Earth is not 
a perfect circle. 

A traditional conception of science regards all physical models as inexact approx- 
imations of reality, and holds that the goal of science is to produce progressively 
more accurate models whose predictions more closely match observations than the 
predictions of previous models. That conception of science is reasonable when the 
values of observable quantities are real numbers, since the real numbers predicted 
by physical models are never exactly the same as the real numbers 'measured' in the 
laboratory. But when non-negative integers are used for the values of observable 
quantities, then an alternate conception of science is possible. 

In this alternate conception there exist faithful models that are in exact agree- 
ment with reality, but perhaps only for a small subset of all physically observable 



^But it should be noted that the model obtained by omitting 02 and T2 from Model [5TT1 is not 
identical to Model [3?T1 In particular, for each state in Model l3Tl there are ten indistinguishable 
states in the model obtained by omitting 02 and T2 from Model 15.11 That is, these models are 
not isomorphic. See Section llOl 
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quantities. For example, Model ISTTl is faithful, but only predicts the angular posi- 
tion of the Earth to within 43.2 degrees, and only for a limited range of times. The 
goal of science is then to produce more refined models. That is, the goal of science 
is to discover faithful models which have larger sets of observable quantities, and 
are therefore capable of predicting increasing numbers of facts. 

6. Radioactive Decay 

Given non-negative integers x and y, let j3{x,y) be the length y sequence com- 
posed of the first y bits in the binary expansion of x. For example /3(13, 6) ~ 
(0, 0, 1, 1, 0, 1). Now suppose that a single atom of a radioactive isotope, such as 
nitrogen-13, is placed inside a detector at time t = Q. We say that the detector has 
status 1 if it has detected the decay of the isotope, and has status otherwise. The 
history of the detector at time t is the length t sequence of bits corresponding to 
the status of the detector at times 1 through t. For example, if the isotope decays 
sometime between t — 2 and < = 3, then the history of the detector at time t — b 
is (0, 0, 1, 1, 1). The following computable physical model models the status of the 
detector as a function of time. 

Model 6.1 (Radioactive Decay). Let S be the set of all triples ^i,/3(2" — l,t),j) 
where n, t, and j are non-negative integers such that n < t, t ^ 0, and 2j < 2" — 1. 
The history of the detector is given by the function ri{t,h,j) = h, and the time, 
measured in units of the half-life of the isotope, is given by the function T{t, h,j) = t. 

This is a model of the many-worlds interpretation |5] of radioactive decay. Sup- 
pose that one asks, "What will the status of the detector be at time t = 2?" There 
are four states {t, h,j) such that r(t, h,i) = 2, namely 

(2,(0,0), 0) (2,(0,1), 0) (2,(1,1), 0) (2,(1,1), 1) 

In three of these states, the detector has status 1, and in one state it has status 
0. If we assume that each state of the system is equally likely, then there is a | 
probability that the detector will have status 1 at time t — 2. But if we ask, "If the 
detector has status 1 at time t — 1, then what will its status be at time t = 2?" The 
answer is "1", since the detector has status 1 at time 2 in both states where the 
detector had status 1 at time 1. These results are in agreement with conventional 
theory. 

7. Ensembles of Physical Models 

Suppose that a planet orbits a distant star and that we are uncertain of the 
planet's orbital period. In particular, suppose that we believe its motion is faithfully 
described by either the discrete model of planetary motion (Model 13. 1|) or by the 
following computable physical model. 

Model 7.1. Let S be the set of all pairs (^{r ; s), (p ; q)) such that 



i 1 I n 1 

p = 360 

10 100 ^ \IQ 100 

i+l I fn+1 1 

^+100 '^ = ''^^+100 



LOoJ) 



n 

10 100- 

1 1 



10 100. 



for some integer i between —20000 and 20000, and such that n = [«/4j • The angular 
position of the planet, represented as a range of angles measured in degrees, is given 
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by the Junction a(^{r ;s),{p;q)') ~ [p'tQ)- The time interval, measured in Earth 
years, is given by the function T(^{r ; s), [p ; q)) ~ {r ]s). 



Note that this model is similar to Model I3.1[ except that the orbital period of the 
planet is 4 Earth years, rather than 1 Earth year. 

If for each of the two models we are given a rational number expressing the 
probability that that model is faithful, then a statistical ensemble of the models 
may be constructed. For example, if Model 13.11 is twice as likely as Model 17.11 
then a corresponding statistical ensemble is the following. Note that this statistical 
ensemble is itself a computable physical model. 

Model 7.2 (Ensemble of Models). Let S be the set of all triples (^{r ;s),{p; q),j) 
such that 



I 

10 



1 



1 
100 

+ 



p = 360 



/ n 1 
llO " 100 ^ 



n 

To 



1 



10 100 
for some integer i between 



q = 360 
-20000 an 



1 



10 



1 
100 



1 

100 

+ 1 



+ 



1 



n 



10 lOOJ 
20000, where j = 0, 1, or 2, and where 

if j — or 1 



LV4J ifj = 2 



The angular position of the planet, represented as a range of angles measured in 
degrees, is given by the function a(^{r ; s) , {p ; q) , j') = (j)\q). The time interval, 
measured in Earth years, is given by the function T(^{r ; s), [p ; q), j) — {r \s). 

Since Model [5TT] is twice as likely as Model lTTTl there are two states, ((r ; s), (p ; g), O) 
and ((r ; s), {p ; q), l) in the ensemble for each state ((r ; s), (p ; g)) in Model [XTl and 
there is one state ((r ; s), (p ; g), 2) in the ensemble for each state ((r ; s), (p ; g)) in 
Model mi Note that the index j in each state {{r ; s), (p ; q), j') of the ensemble is 
not observable. 

Now, if we ask for the position of the planet during the time interval (0.29 ; 0.41), 
for example, there are three possible states ((r ; s), (p ; g), j) in the ensemble such 
that 

T((r;s),(p;g),j)- (0.29; 0.41) 

namely 

((0.29; 0.41), (104.4; 147.6), 0) 

((0.29; 0.41), (104.4; 147.6), l) 

((0.29; 0.41), (356.4; 39.6), 2) 

Since the planet's angular position is (104.4 ; 147.6) for two of these three states, 
the position measurement (104.4; 147.6) has a probability of |. Similarly, because 
the planet's angular position is (356.4 ; 39.6) for one of the three states, the po- 
sition measurement (356.4; 39.6) has a probability of ^. These probabilities are 
a direct reflection of our uncertainty about which of the two underlying physical 
models. Model \3A\ or Model [7?T1 is the true faithful model. In particular, because 
Model l3TT] has been deemed twice as likely as Model FTTl the position of the planet in 
Model IXn namely (104.4 ; 147.6), has twice the probability of the position predicted 
by Model [m namely (356.4 ; 39.6). 
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It is important to note that there is no observable quantity corresponding to 
probabihty in Model 17.21 Instead, probability is a mathematical tool used to inter- 
pret the model's predictions. This sort of interpretation of an ensemble of models 
is appropriate whenever the ensemble is composed from all possible models which 
could describe a particular system, with the number of copies of states of the in- 
dividual models reflecting our confidence in the predictions of those models. See 
reference |10j for a more detailed account of this subjectivist interpretation of prob- 
ability in physics. 

Ensembles may be constructed in other circumstances as well, and we may refer 
to such ensembles as non- statistical ensembles of physical models. Non-statistical 
ensembles of physical models are commonplace in the sciences. For example, they 
result whenever a constant, such as an initial position, is left unspecified in the 
statement of a model. That model can then be used to describe any member of a 
family of systems, each of which may have a different value for the constant. But 
most importantly, when a non-statistical ensemble of physical models is constructed, 
no claims as to the likelihood of one value of the constant, as compared to some other 
value of the constant, are being made. In fact, this is the defining characteristic of a 
non-statistical ensemble of models. Non-statistical ensembles can be useful because 
they provide a convenient way to collect together sets of closely-related models. 

8. Incompatible Measurements 

A pair of measurements is said to be simultaneous if and only if they are both 
performed while the system is in a single state. An essential feature of quantum 
mechanical systems is that there may be quantities which are not simultaneously 
measurable. For example, the measurement of one quantity, such as the position 
of a particle, might affect the subsequent measurement of another quantity, such 
as the particle's momentum. Such measurements are said to be incompatible. It is 
natural to ask whether computable physical models can be used to describe systems 
which feature incompatible measurements. 

Discrete quantum mechanical systems are often formalized as follows [6l[24]. The 
quantum mechanical state of a system is a normalized vector v in some normed 
complex vector space V. Typically, V is a Hilbert space and is a wave func- 
tion. For each quantum mechanical measurement there is a corresponding set 
B = {vi,V2,V3, . . .} of normalized basis vectors for V. Each member of B cor- 
responds to a possible value of the measurement. Because i3 is a basis for V, 
V = aivi + a2V2 + a^vs + ■ ■ ■ for some complex numbers ai, 02, 03, .... If the 
system is in quantum mechanical state v and no two members of B correspond to 
the same measurement valueQ then the probability that the measurement will have 
the value corresponding to w„ is la^p. In this case, if the actual value which is 
measured is the value corresponding to u„, then the quantum mechanical state of 
the system immediately after that measurement is Vn. The state v is said to have 
collapsed to w„. During the time between measurements, the quantum mechanical 
state of a system may evolve according to a rule such as Schrodinger's equation. 



^Alternatively, if Vm , f ^2 i > ■ • • ^re distinct basis vectors corresponding to the same mea- 
surement value, then the probability of measuring the value is \ani |^ -I- \a,n2 |^ -I- |i>i3 + ■ ■ ■ ■ If 
that value is actually measured, then the state of the system immediately after the measurement 
is the normalization of a„j^v„-^ + o.„2V„2 + '^n^Vn^ See reference |13| . 
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Consider, for example, the problem of measuring the components of the spin 
of an isolated electron. In this case, V is the set of all vectors (a, b) such that a 
and b are complex numbers, where the norm ||(a, is defined to be y^aP + |6p. 
A quantum mechanical measurement of the z component of the electron's spin 
has two possible values, —^fi and +^fi. The basis vectors corresponding to these 
values are (0,1) and (1,0), respectively. The quantum mechanical measurement 
of another component of the electron's spin, lying in the xz plane at an angle of 
60 degrees to the z axis, also has two possible values, —5 ft and +-^h. The basis 

vectors corresponding to these values are (— 5, -^) and (-^, respectively. So, 
for example, if the spin component in the z direction is measured to have a value 
of at time t = 0, then since 

(>.»)-K4-f)-f(#^) 

there is a |-^| = | probability that if the 60-degree electron spin component is 
measured at time t = 1, then that component will also have a value of +-^h. 

Supposing that the 60-degree electron spin component is measured to have a 
value of +-^h at time i = 1, a similar line of reasoning implies that if the spin's z 
component is measured at time t — 2, then there is a — j probability that the 
value of that measurement will be — ^ft, since 

(f4)4<°'')-f<^'°) 

Therefore, if the z component of the electron's spin is measured at time i = 0, 
followed by a measurement of the 60-degree spin component at time t = 1, and 
followed by another measurement the z component at time t = 2, then the values 
of the two measurements of the z component need not be the same. Indeed, the 
quantum mechanical state of the system does not chang^ between times t = and 
t — 1, 01 between times t = 1 and t = 2, but the measurement of the 60-degree 
spin component at time t — 1 disturbs the system and can potentially change the 
value of any subsequent measurement of the z component. That is, measurement 
of the electron's 60-degree spin component is incompatible with measurement of its 
z component. 

Let us formalize this system as a computable physical model. The system is 
composed of the electron, the apparatus used to make the quantum mechanical 
measurements, and the researcher who chooses which components to measure We 
assume that the quantum mechanical state of the electron is (1,0) at time t = 0, 
and that the researcher makes subsequent quantum mechanical measurements of the 
electron's spin components at times t — 1 and t — 2. When a quantum mechanical 
measurement is performed, a record is made (perhaps in the researcher's notebook) 
of the value of this measurement and of the component that was measured. We 



°In this case, the quantum mechanical state of the system does not change between measure- 
ments because the electron is isolated. For example, the electron is free from external electromag- 
netic fields or other influences that might cause its spin to precess. 

'^We refrain from asking questions about the probability with which the researcher chooses 
which components to measure. That is, this model describes a non-statistical ensemble of 
researchers. 
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construct the computable physical model of this system from the point of view of 
an agent who observes only this recorded history and the time. 

Model 8.1 (Electron Spin Measurement). Let S be the set of all triples {t,h,j) 
such that 

t=l t=l t=l 

h={0;+l) or /i = (60;-l) or /i = (60;+l) 

j = j = 1 j = 2 + m 



or 



t=2 t=2 

/i = ((0;+l),(0;+l)) or ft = ((0 ; +1), (60 ; -1)> 

j = 5 j = 6 



or 



t = 2 

/i = ((0;+l),(60;+l)) 
j = 7 + m 



t = 2 

/i=((60;-l),(0;-l)) 
j = 10 



or 



t = 2 

/i = ((60;-l),(0;+l)) 
7 = 11 + m 



t = 2 

or /i = ((60;-l),(60;-l)> 

j = U 



or 



t = 2 t 
ft = ((60;+l),(0;-l)) or h 

j = 15 + n j 

or 

t = 2 

/i = ((60;+l),(60;+l)) 
i = 27 + m 

for some integers m and n with < m < 2 and < n < 8. The time is given 

by the function T{t.h, j) = t. The history is given by the function rj{t^h,j) — h. 
A history is a chronological sequence of records, with the leftmost record being the 
oldest. Each record is a pair {a ; b) of rational numbers, where a is the angle from 
the z axis, measured in degrees, of a component of the electron's spin, and where b 
is the value of that component, measured in units of ^H. 

Note that each state {t,h,j) has a distinct index j, which we will use to identify 
that particular state. 



= ((60;+l),(0;+l)) 
= 24 + rn 
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Model \8A] corresponds to the quantum mechanical system in the following sense. 
First, the quantum mechanical state of the system at time t corresponds to a set 
of states in the computable physical model. For example, if the researcher decides 
to measure the 60-degree component of the electron's spin at time t = 1, then the 
quantum mechanical state of the system is represented by the set of states with 
indices 1 through 4. Assuming that the states in the set are equally likely, there 
is a I probability that this component will have a measured value of +5^1, for 
example. Immediately after the measurement is made, the quantum mechanical 
state collapses, becoming either the set of states with indices 2 through 4, or the 
singleton set containing only the state with index 1. The collapse occurs because the 
information provided by the quantum mechanical measurement allows us to identify 
the state of the system more precisely, eliminating those states which disagree with 
the measurement result 13 The quantum mechanical state then evolves to a new set 
of states at time t = 2. For example, if the measured value of the 60-degree electron 
spin component is at time t = 1, and if the researcher plans to measure the 
0-degree electron spin component (that is, the z component) at time t = 2, then 
the quantum mechanical state immediately before that measurement at time t — 2 
is the set of states with indices 15 through 26. 

Computable physical models similar to Model 18.11 can be constructed for quan- 
tum mechanical systems which satisfy the following criteria. 

(1) There is a set of possible measurements {mo, mi, m2, . . . , m^, . . .} indexed 
by non- negative integers i. 

(2) Every discrete time step, one measurement from this set is performed. 

(3) The possible values of each measurement m,; are identified with non- negative 
integers. 

(4) If (f){i, n, t, h) is the probability that the measurement with index i has the 
value n, given that the measurement is performed at time step t and that 
h = ((ii ; ni), . . . , ; nt_i)) is the history of past measurements and 
their values, then (/>(i, n, t, h) is a rational number. 

(5) If there is no measurement with index i or if the non-negative integer n 
does not correspond to a value of the measurement with index i, then 
4>(i, n, t, h) = 0. 

(6) For each choice of non-negative integers i, t, and h, there are only finitely 
many non-negative integers n such that 4>{i, n, t, h) > 0. 

(7) (/) is a recursive function. 

If a quantum mechanical system satisfies these criteria, then we can determine 
whether or not 

s = {t,{{ii;ni), (i2 ;n2), . . . , {it;nt)),j) 
is in the set S of states of the corresponding computable physical model as follows. 
First, if t = 0, then s is not in S. Next, let hi — and for each positive integer k 
with 1 < k <t,let 

hk = {{ii ; ni), (12 ; 71-2), • ■ • , («fc-i ; »^fe-i)) 
Now we perform the following calculations for each positive integer k < t. If 
</>(*fc, n-k, k, hk) — 0, then s is not in S. Otherwise, there must be finitely many non- 
negative integers n such that the probability (j)(ik,n,k,hk) is greater than zero. 



For a more detailed discussion of this ensemble interpretation of the collapse of a quantum 
mechanical state, see reference [l]. 
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Since probabilities must sum to 1, those values for n may be found exhaustively by 
calculating (j){ik, 0, fc, hk)-, (piik, 1, k, hk), (f>{ik, 2, fc, hk), and so on, until the the sum 
of these probabilities reaches 1. Let dk be the least common denominator of these 
rational probabilities, and let be the unique positive integer such that 

(t){ik,nk,k,hk) = ^ 
dk 

If j < aia2 ■ ■ ■ at, then s is in S. Otherwise, s is not in S. 

9. The Accuracy of Measuring Instruments 

An important feature of the discrete model of planetary motion (Model 13. ip is 
that the intervals representing time and angle measurements overlap. The amount 
of overlap between adjacent intervals is determined by the accuracy of the corre- 
sponding measuring instrument. The introduction of overlapping intervals is moti- 
vated by an argument such as the following. 

If Model [??T] were constructed using disjoint, non-overlapping intervals, then the 
states ((r ; s), (p ; g)) of that model would be given by 

r^i/lQ p = 360(r-[rj) 

s = (i-|-l)/10 g = 360(s-[sj) 

where i is an integer. In particular, ((0.2 ; 0.3), (72 ; 108)) and ((0.3 ; 0.4), (108 ; 144)) 
would be two such states, with (r ; s) representing the state's time interval, mea- 
sured in years, and with (p ; q) representing the corresponding interval of angular 
positions for the Earth, measured in degrees. According to this model, if the posi- 
tion of the Earth is measured at time t = 0.298 years, then t is within the interval 
(0.2 ; 0.3), and the state of the system is ((0.2 ; 0.3), (72 ; 108)). Therefore, accord- 
ing to this model, the position of the Earth should be between 72 and 108 degrees. 
Indeed, the simple model of planetary motion (Model [TTT|) predicts that the angular 
position of the Earth at time t = 0.298 years should be 360(0.298 - [0.298J ) « 107 
degrees. But the true position of the Earth in its orbit deviates from Model 11.11 
In this case, the true position of the Earth at time t = 0.298 years is about 109 
degrees which is outside the interval (72 ; 108). Therefore, if the discrete model 
were constructed using disjoint, non-overlapping intervals, then the model would 
fail when t = 0.298 years. 

But the discrete model of planetary motion (Model 13. 1|) was constructed using 
overlapping intervals. In particular, 

((0.19 ; 0.31), (68.4 ; 111.6)) ((0.29 ; 0.41), (104.4 ; 147.6)) 

are two states in Model [01 Note that at time t — 0.298 years. Model [3Tl could be 
in either of these two states. Furthermore, any pair of real-valued time t and angle 
a measurements which satisfy 

\a-im{t~ \ t\) \ < 7.2 

fall within the time and angle intervals of some common state of Model 13.11 Since 
\a—360{t— [t\)\ is at most 2 degrees [7] for all physically observed angles a measured 
at times t, Model [3Tl is faithful. 

-'^-'^This is assuming that time is measured in anomalistic years, with each year beginning at 
perihelion passage. During the course of a year, the position of the Earth is the true anomaly, 
measured relative to that perihelion passage. 
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It is important to point out, though, that Model [Ql is faithful only if the results 
of measurements are uncertain when they occur within the region of overlap. For 
example, at time t = 0.306 years, two results of a time measurement are possible, 
(0.19; 0.31) and (0.29; 0.41), and an observer cannot be certain which of these 
intervals is the value of the measurement. The actual angular position of the Earth 
at time t = 0.306 years is about 112 degrees, so (104.4; 147.6) is the only possible 
result of a position measurement. Since 

((0.19; 0.31), (104.4; 147.6)) 

is not one of the states of Model l3.11 the observer is expected to realize, in retrospect, 
after measuring the angular position, that the true time measurement must have 
been (0.29 ; 0.41). After providing a model for the phenomenon of accuracy, we will 
be able to reformulate Model 13.11 so that the results of measurements no longer 
possess this sort of ambiguity. 

But first, note that the accuracy of a measuring instrument, by definition, can 
only be quantified relative to some other, more precise quantity. For example, 
the argument above, concerning accuracy in Model 13.11 makes frequent reference 
to exact real-valued angles and times. Indeed, even when we express an angle 
measurement as an interval, such as (68.4 ; 111.6), we are implying that it is possible 
to distinguish an angle of 68.4 degrees from an angle of 111.6 degrees, and that other 
angles lie between those two values. In principle, though, it is possible to describe 
the accuracy of a measuring instrument in a purely discrete manner, without any 
mention of real numbers. For example, let us consider an instrument for measuring 
distances in meters, with the value of a measurement represented as an integer 
number of meters. The accuracy of this measuring instrument can be quantified 
relative to a second instrument which measures distances in decimeters. 

Presumably, the phenomenon of accuracy results from our inability to properly 
calibrate measuring instruments. Although there are many different underlying 
causes of calibration error, it suffices to consider only one such cause for a simple 
model of this phenomenon. We will suppose that when we measure a distance in 
meters, that we have difficulty aligning the measuring instrument with the origin, 
so that sometimes the instrument is aligned a decimeter too far in the negative 
direction, and at other times a decimeter too far in the positive direction. Hence, 
there are two different physical models for the measurement. In one model the 
instrument is misaligned in the negative direction, and in the other model it is 
misaligned in the positive direction. Since we do not know which of these two 
models describes any one particular measurement, it is appropriate to combine 
them in the following statistical ensemble. 

Model 9.1. Let S be the set of all triples (^(m), C((i), C(i)) such that 

d + i\ 



where d is an integer, and where i = —1 or +1. The distance, measured in meters, 
is given by the function 

/i(C(m),C(d),C(0) =C(m) 
The same distance, measured in decimeters, is given by the function 

5(C(m),C(d),C(0) -C(rf) 
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Note that the function ( was defined in Section [2] Also note that the index i 
in each state (C(m), represents the cahbration error, which is cither — 1 

decimeter or +1 decimeter. Model IHTTl is a computable physical model. 

A measurement of d decimeters in Model \9A\ can be interpreted as correspond- 
ing to an interval of ; ^j^) meters. Note that a measurement of 9 decimeters 
(corresponding to an interval of (0.9; 1.0) meters), for example, is possible in two 
distinct states of the model: 



(C(0),C(9),C(-i)) 



(C(i),C(9),C(+i)) 



Similarly, a measurement of 10 decimeters (corresponding to an interval of (1.0 ; 1.1) 
meters) is possible in the states 



(C(o),C(io),C(-i)) 



(C(i),C(io),C(+i)) 



Hence, a measurement of meters overlaps with a measurement of 1 meter on 
the intervals (0.9 ; 1.0) and (1.0 ; 1.1). And in general, a measurement of m meters 
overlaps with a measurement of m + 1 meters on the intervals (m + 0.9 ; m + 1.0) 
and (m + 1.0 ; m + 1.1). Therefore, a measurement of m meters in Model \9A\ can 
be understood as corresponding to an interval of (m — 0.1 ; m + 1.1) meters, with 
adjacent intervals overlapping by 0.2 meters. 

Of course, this interpretation of Model 19.11 presumes that decimeters can be 
measured with perfect accuracy. A more realistic computable physical model can be 
constructed by supposing that decimeter measurements can also be misaligned, for 
example, by —1 centimeter or +1 centimeter. Note that centimeters are treated as 
unobserved, purely theoretical constructions in this model — there is no observable 
quantity for centimeter measurements. 

Model 9.2. Let S be the set of all quintuples (C(m), C('^); C(c), C(*)i C(j)) such that 



lOi 



100 



d = 
-I or +1 



c + j \ 
10 J 

. The distance, measured 



where c is an integer, i = —1 or +1, and j 
in meters, is given by the function 

M(C(m),C(d),C(c),C(0,C(j)) =C(m) 

The same distance, measured in decimeters, is given by the function 

S{C{m), C{d), Cic), m, Cij)) = (id) 

As an important application, the model of accuracy described in this section can 
be used to reformulate the discrete model of planetary motion (Model [O]) . 



Model 9.3. Let S be the set of all quintuples ((r ; s), (p ; g), C(i), C(j'), ^(fc)) such 
that 

k + j\ 



k + i 
10 



m 

To 



10 



1 

100 

1 1 

~ ^ 100 



10 J 
360 

360'"" 



1 
100 



1 



10 



1 

100 



1 

^ 100 
n + 1 
10 



1 

100. 
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for some integers m and n, for some integer k between —200000 and 200000, and 
where i = —1 or +1, and j = —1 or +1. The angular position of the Earth, 
represented as a range of angles measured in degrees, is given by the function 

a{ir■s),ip■,q),Ci^),aJ)Xik))^ip■,q) 

The time interval, measured in years, is given by the function 

r((r;s),(p;q),C«,CO-),C(fc)) = (r;s) 

Note that like Model lO] this model is faithful. But the faithfulness, in this case, no 
longer requires that the results of some measurements be uncertain. Instead, given 
any particular measurement, the state of the system is uncertain. For example, 
there are 40 distinct states s such that r(s) = (0.29 ; 0.41). 

In contrast to Model 13.11 consider what happens if Model 19.31 is used to explain 
measurements taken at time t — 0.306 years. Associated with each state 

((r;s),(p;g),CW,C(j),C(fc)) 

in Model l973l is an integer k, intended to represent the time interval (j^q ' ToS") 
during which the system is in that state. At time t = 0.306 years, k — 30, and the 
system could be in one of the following four states: 

((0.19 ;0.31), (68.4; 111.6), C(-1),C(-1),C(30)) 
((0.29 ; 0.41), (68.4 ; 111.6), C(+l), C(-l), C(30)) 
((0.19 ; 0.31), (104.4 ; 147.6), C(-l), C(+l), C(30)) 
((0.29 ; 0.41), (104.4 ; 147.6), C(+l), C(+l), C(30)) 

Like Model [331 two time measurements are possible, (0.19; 0.31) or (0.29; 0.41). 
And since the actual angular position of the Earth at time t = 0.306 years is about 
112 degrees, the measured position of the Earth is (104.4; 147.6) degrees at that 
time. Unlike Model I3TT1 this position measurement is compatible with either time 
measurement, since 

((0.19 ; 0.31), (104.4 ; 147.6), C(-l), C(+l), C(30)) 

((0.29 ; 0.41), (104.4 ; 147.6), C(+l), C(+l), C(30)) 

are both states of Model 19.31 



10. Isomorphism Theorems 

Given a physical model with a set S of states and a set A — {ai, a2,a^, . . .} of 
observable quantities, we write {S, A) as an abbreviation for that modelri 

Definition 10.1. Two physical models (5*,^) and {T,B) are isomorphic if and 
only if there exist bisections cf) : S T and ijj : A ^ B such that a{s) — tpio') {'Pis)) 
for all s € S and all a d A. 



In the interest of generality, the definition of a computable physical model places no restric- 
tions on the set A except that its members must be total recursive functions. Some authors prefer 
to restrict their attention to finite sets A. For example, see reference [3]. Other authors may 
prefer to restrict their attention to observable quantities computed by programs which belong to 
a recursively enumerable set. 
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Intuitively, isomorphic physical models can be thought of as providing identical 
descriptions of the same system^ 

Given any particular computable physical model {S, A), there are many different 
models which are isomorphic to (S^A). The following two theorems provide some 
convenient forms for the representation of computable physical models. Let tt" be 
the projection function that takes a length n sequence of non-negative integers and 
outputs the ith element of the sequence. That is, Trf (xi, 2:2, . . . , x„) = Xi for any 
positive integer i < n. 

Theorem 10.2. If A is a finite set, then the computable physical model (S, A) is 
isomorphic to some computable physical model whose observable quantities are all 
projection functions. 

Proof. Given a computable physical model (S^A) with A = {ai,Q;2, . . . ,Q;„}, let 
(T, B) be the computable physical model such that 

T = { (ai(s),a2(s), . . . ,a„(s),s) | s G S" } 

and let B be the set of projection functions {tt"'*'^, ttj"*^^, . . . ,7rJJ+^}. By construc- 
tion, {T,B) is a computable physical model isomorphic to {S,A). □ 

Theorem 10.3. If S is an infinite set, then the computable physical model iS,A) 
is isomorphic to some computable physical model whose set of states is the set 
of all non-negative integers. If S has n elements, then the computable physical 
model (S, A) is isomorphic to some computable physical model whose set of states 
is {0, 1, . . . , n — 1}. 

Proof. By definition, if {S, A) is a computable physical model, then 5 is a recursive 
set. It immediately follows that S is recursively enumerable. In particular, if S is 
infinite, then let T be the set of non-negative integers and there is a bijective recur- 
sive function ip from T to S. If S has n elements, then there is a bijective recursive 
function ij; from T = {0, 1, . . . , n — 1} to S. Now, given A = {ai, a2, as, . . .}, let 
B = {aioip, a20'4>, a^oip, . . .}, where aoip denotes the composition of the functions 
a and -0- By construction, {S,A) is isomorphic to {T,B). □ 

Although Theorem 110.31 implies that any computable physical model (5*, A) is 
isomorphic to a computable physical model (T, B) where T is a set of consecu- 
tive non-negative integers beginning with zero, there is no effective procedure for 
constructing a program that computes the characteristic function of T, given a pro- 
gram for computing the characteristic function of S when 5* is a finite set. That is. 
Theorem 110.31 does not hold uniformly. 

Definition 10.4. A non-negative integer physical model is a pair {S, A) where S 
is a set of non-negative integers and where each member of A is a partial function 
from the non-negative integers to the non-negative integers. S is the set of states 
of the model, and A is the set of observable quantities of the model. 

Note that if a is an observable quantity of a non-negative integer physical model 
{S,A), then a might be undefined for some inputs. The non-negative integer phys- 
ical models form a more general class of objects than the computable physical 



^•^Rosen 1191 defined a weaker notion of isomorphism. Physical models that are isomorphic in 
Rosen's sense are not necessarily isomorphic in the sense described here. 
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models. In particular, a computable physical model is a non-negative integer phys- 
ical model whose set of states is a recursive set and whose observable quantities are 
total recursive functions. 

Theorem 10.5. A non-negative integer physical model {S, A) is isomorphic to 
some computable physical model if S is a recursively enumerable set and if each 
member of A is a partial recursive function whose domain includes all the members 
ofS. 

Proof. Given a non-negative integer physical model {S, A) , note that the construc- 
tion of the computable physical model (T, B) in the proof of Theorem 110.31 only 
requires that 5* be a recursively enumerable set and that each member of A be a par- 
tial recursive function whose domain includes all the members of S. Therefore, any 
such non-negative integer physical model (S*, A) is isomorphic to the computable 
physical model (T, _B). □ 

Let (f) be any partial recursive function and let S be the largest set of consecutive 
non-negative integers beginning with zero such that (j){s) is defined for each s £ S. 
Let A be the set {tt™ o 0, tt™ o 0, . . . , tt™ o ((>} for some non-negative integer m. By 
Theorem 110.51 {S,A) is isomorphic to a computable physical model. We say that 
any such computable physical model is determined by (j). 

Theorem 10.6. Every computable physical model with finitely many observable 
quantities is determined by some partial recursive function (j). 

Proof. Let {S, A) be a computable physical model with A — {ai, a2, . . . , am}- Let 
ip be the recursive function given in the proof of Theorem 110.31 If S has only n 
states, then let </>(«) be undefined for all non-negative integers i > n. Otherwise, 
define 

(P{i) = (ai('0(j)), a2{ip{i)), a„i{il;{i))) 
By construction, (S, A) is determined by 0. □ 

A physical model {S, A) is said to be reduced if and only if for each pair of distinct 
states si and S2 in 5, there exists an a € A with a{si) ^ a{s2). 

Theorem 10.7. If{S,A) and (T,B) are isomorphic physical models and {S,A) is 
reduced, then (T, B) is also a reduced physical model. 

Proof. Let (S*, A) and (T, B) be isomorphic physical models and let (5, A) be re- 
duced. Since {S, A) and (T, B) are isomorphic, there exist bijections (f> : S ^ T and 
tjj : A^ B such that a{s) = ^(a)(0(s)) for all s £ 5 and all a G A. Now suppose 
that ti and ^2 are distinct states in T. Because ^ is a bijection, (j)~^{ti) and (/)^^(i2) 
are distinct states in S. But S is reduced, so there must exist an a G A such that 
a((/>-i(ii)) a{(l)-^{t2)). Furthermore, 

and 

a{r\t2)) = i;{a){4>{r\t2))) = ^(a)(t2) 

Hence, there exists a f3 & B such that (3{ti) ^ /3(t2), namely /3 = We may 

conclude that the physical model (T, B) is reduced. □ 
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An epimorphism from a physical model (S", A) to a physical model (T, B) is a 
pair of functions (</>, -0) such that is a surjection from S to T and is a bijection 
from A to B, where a{s) = ^{a)[(f){s)) for all s e S* and all a £ A. Two physical 
models {Si, Ai) and (52, A2) are said to be observationally equivalent if and only if 
there are epimorphisms from {Si, Ai) to (T, B) and from (5*2, A2) to {T, B), where 
(T, B) is some reduced physical model. 

Theorem 10.8. // {Si,Ai) and (52,^2) are isomorphic physical models, then 
{Si,Ai) and (52,^2) are observationally equivalent. 

Proof. Define an equivalence relation on ^2 so that r G 6*2 is related to s G 6*2 if and 
only if a{r) — a{s) for all a G A2. Let T be the corresponding set of equivalence 
classes of 52- For each a ^ A2, define a function a' so that ii s G t G T, then 
a'{t) = a{s). Let B = {a' \ a ^ A2}. By construction, {T,B) is a reduced 
physical model. Also note that there is an epimorphism {(l),ip) from (52,^2) to 
(T, B). Namely, is the function that maps each member of ^2 to its corresponding 
equivalence class in T, and ip is the function that maps each a d A2 to a' B. 

Now suppose that (5*1,711) and (52,^2) are isomorphic. By definition, there are 
bijections (f>' from 5i to 52 and ip' from Ai to A2 such that a{s) = ■0'(a) ((/>'(s)) for 
all s € 5i and all a (z Ai. Since (0, ip) is an epimorphism from (52, A2) to (T, B), it 
immediately follows that (0 o 0', ?/) o ?/;') is an epimorphism from (5i, Ai) to (T, B). 
We may conclude, by definition, that (5i,Ai) and (52,^2) are observationally 
equivalent. □ 

Intuitively, two physical models are observationally equivalent when they both 
make the same observable predictions. For example, as was discussed in Section [SJ 
the model obtained by omitting the observable quantities a2 and T2 from Model [SHI 
is observationally equivalent to the discrete model of planetary motion (Model [XT]) . 
Moreover, if a physical model {Si,Ai) is faithful, and if (5i, Ai) is observationally 
equivalent to (52,^2), then (52,^2) is also faithful. 

It is important to note that the converse of Theorem 110.81 does not hold. That 
is, observationally equivalent models are not necessarily isomorphic. Consider, for 
example, the computable physical models ({0,1}, {a}) and ({0, 1, 2}, {/?}) where 
a(s) = s for ah s e {0,1} and where /3(s) = [s/2j for ah s G {0,1,2}. These 
models are not isomorphic because {0, 1} and {0, 1, 2} have different cardinalities. 
Yet, they are observationally equivalent, since both models have a single observable 
quantity whose only possible values are and 1. Physical models that are isomor- 
phic must not only make the same observable predictions, but must also have the 
same structure. The models ({0,1}, {a}) and ({0, 1, 2}, {/?}) have different struc- 
tures because, assuming that the states are equally likely, they both give different 
answers to the question, "What is the probability that the observable quantity has 
value 0?" 

11. Oracles and Effective Topologies 

Given any set X, we can impose a topology on X. Let B he a. basis for this 
topology. The members of B are said to be basis elements. We say that a set 
£3; C S is a local basis for a point x Cz X it and only if the following two conditions 
hold. 

(1) For each L ^ £x, x is a, member of L. 

(2) For each B G B with x G B, there exists a.n L E Cx with L C B. 
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Note that every point x ^ X has a local basis. For example, the set 

C^^ {B C,B\x e B} 

of all basis elements that contain a; is a local basis for x. 

If the basis B is countable, then each basis element can be encoded as a non- 
negative integer. In that case, choose some encoding and let v{n) be the basis 
element encoded by n. We allow for the possibility that a basis element may be 
encoded by more than one non-negative integer. (That is, v is not necessarily an 
injection.) The domain of v, denoted domg v, is the set of all non-negative integers 
n such that v{n) G B. For any function (f) : A ^ B and any set CCA, let 
— {0(c) I c G C} denote the image of C under (f>. We let N denote the set of 
non-negative integers. 

Definition 11.1. A function : N — > domg ^' is said to be an oracle for a point 
X, with basis B and coding v, if and only if z^((/)(N)) is a local basis for x. 

An oracle 4> for x is complete if and only if every n £ domj? v such that x G ^{n) is a 
member of </'(N). An oracle is said to be nested if and only if v{(j){n+V)) C i/{(j)[n)) 
for all n G N. 

A pair (B, v) is said to be an effective topology if and only if S is a countable basis 
for a To topology and v is a, coding for B. Effective topologies were first introduced 
in the theory of type-two effectivity [121 US] . In accordance with that theory, we 
use an oracle for x, with a basis B and coding v, as a representation of the point x 
in an effective topology {B, v). Because effective topologies are Tq, no two distinct 
points are ever represented by the same oracle. 

Of special interest are effective topologies where the subset relation 

{ {bi, 62) I Q v{b2) & 61 G dome & 62 G domg v } 

is a recursively enumerable set0 In particular, if (S, v) has a recursively enumer- 
able subset relation, then 

dome v — {b \ v{b) C v{h) & 6 G dome } 

is also a recursively enumerable set. 

Theorem 11.2. Let (p be an oracle for x in an effective topology {B, i^) with a 
recursively enumerable subset relation. Then there exists a complete oracle ip for x 
in {B, v) that is recursive relative to (j) uniformly. 

Proof. Suppose that cj) is an oracle for x in an effective topology {B, v) with a 
recursively enumerable subset relation. Since z^((/)(N)) is a local basis for x, it 
follows that for each b G domg x & v{b) if and only if there exists an n G N such 
that i^((?I>(n))c v[b). Hence, 

{ 6 G dome | (3 71 G N) [^{(lyin)) C v{b)\ ] 

is the set of encodings of all basis elements that contain x. But this set is recursively 
enumerable relative to because {B, v) has a recursively enumerable subset relation. 
Therefore, there exists a function i/j : N — > domg v, recursive relative to (j), such 
that i/(-!/'(N)) is this set. By definition, t/j is a complete oracle for x. □ 

effective topology with a recursively enumerable subset relation is an example of a com- 
putable topology, as defined in reference I26| . Not all computable topologies have recursively 
enumerable subset relations. 
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Theorem 11.3. Let 4> be an oracle for x in an effective topology {B,v) with a 
recursively enumerable subset relation. Then there exists a nested oracle ip for x in 
v) that is recursive relative to (j) uniformly. 

Proof. Suppose that is an oracle for x in an effective topology (i3, v) with a 
recursively enumerable subset relation. Note that for each pair of basis elements Bi 
and B2 such that x G Bir\B2, there exists a basis element B3 with x G i?3 C BiC\B2, 
by the definition of a basis. Therefore, since z^(0(N)) is a local basis for x, there 
must exist an m G N such that 

X e v{4){m)) c B3 c Bi n B2 

Now define -0 : N — dome v recursively, relative to 0, as follows. Let 0(0) = 0(0) 
and for each n g N let 0(n + 1) = 0(to) for some m G N such that 

v[(j){m)) C v[ip{n)) n v{(t){n + 1)) 

We can find m recursively given 0(n) and because the subset relation for [B, v) 
is recursively enumerable, and the set of all m G N such that 

v{4>{m)) C i/(0(ri)) & v[(t>{m)) C zy(0(n + 1)) 

is therefore recursively enumerable relative to 0. We may conclude that ip is nested 
because 

v^il^in + 1)) C v{i:{n)) n i^(0(n + 1)) C y{ilj{n)) 
for all n e N, and that ^ is an oracle for x because x £ i^{^ip{0)) = i^(0(O)) and 

X € v (0^(n + 1)) C v{i^{n)) n v{(l){n + 1)) C v((j){n + 1)) 
for all n e N. □ 

Define t(a ; &) to be the set of all real numbers x such that a < x < b, and 
let X be the set of all t(a ; b) such that a and b are rational numbers with a < b. 
The members of I are said to be rational intervals. Note that I is a basis for 
the standard topology of the real numbers. Indeed, the oracles for real numbers 
that were introduced in Section |4] were nested oracles with basis I and coding t. 
Another basis for the standard topology of the real numbers is the set Iio,c of 
decimal intervals with accuracy factor c, where c is a positive rational number, and 
where Iio,c is defined to be the set of all t(a ; b) such that 

m c ^ m + 1 c 

"~I0"~10" " 10" ^ 10" 

for some integer m and some positive integer n. We call n the number of digits of 
precision of {a;b). 

Note that both (Z, t) and (Iio,c , '-) have recursively enumerable subset relations. 
The following theorem asserts that if is an oracle for a real number x with basis 
X and coding t, then there exists an oracle ip for x with basis Xio,c and coding t 
that is recursive relative to uniformly. 

Theorem 11.4. Let {Ajv) and {B,v) be effective topologies with recursively enu- 
merable subset relations such that B Q A, and such that A and B are bases for the 
same topology. If is an oracle for x in {A^ v), then there exists an oracle "0 for x 
in (B, i>) that is recursive relative to uniformly. 
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Proof. Suppose that {A, v) and (S, v) are effective topologies as described in the 
statement of the theorem, and that (p is an oracle for a point x in [Atv). By 
definition, 

for every n G N. And because i^(0(n)) is an open set, it is a union of basis elements 
from B. Hence, there must exist a B E B such that 

X e B <Z v{^(j}{n)) 

But J/(</)(N)) is a local basis for x, and _B is a basis element in A, so there exists an 
TO G N such that x £ v {(j){m)) C B. Therefore, we have that for each n € N there 
exist B E B and m G N such that 

X G v{(j){m)) CSC v{(j}{n)) 

Now, since (B, v) has a recursively enumerable subset relation, domg v is & re- 
cursively enumerable set. Then, because (.4, v) also has a recursively enumerable 
subset relation, the set 

{ (m, b) I z^((/)(to)) C v{h) C v{(j){n)) &i m £ N k b e dome } 

is recursively enumerable relative to 0, for any n G N. Therefore, there is a function 
■0 : N — > dome f^, recursive relative to 0, such that ^(n) — b for all n G N, where 

j/((/)(to)) C iy{b) C 

for some m G N. But this function ^ is an oracle for x in (S, t^), because 

X G v{(t){rn)) C i^(?/;(n)) C :/((/)(n)) 
for all 71 G N. □ 

12. Basic Representations of Sets 

Definition 12.1. Let {B^v) be an effective topology on a set X, and let A be any 

subset of X . We say that a set R of non-negative integers is a basic representation 
of A in the effective topology {B^v) if and only if the following two conditions hold. 

(1) R C dome v 

(2) X £ A if and only if there exists a local basis for x with C u[R). 

Note that condition [2] of the definition ensures that no two distinct sets in (S, v) 
have the same basic representation. Note further that if i? is a basic representation 
of A, then {AC\ v{r) | r G i? } is a basis for the subspace topology on A, and this 
is an effective topology with coding \r\_A n v{r)\ . 

In an effective topology we use basic representations to represent sets of points, 
but not all sets of points have basic representations. For example, there are 2^ " 
many sets of real numbers, but since a basic representation is a set of non-negative 
integers, there are at most 2^'' many basic representations. Nevertheless, many 
commonly-studied sets have basic representations 

"'^^In the effective topology (X, t), the set of rational numbers does not have a basic represen- 
tation, but the set of irrational numbers has the basic representation 

I. Vn! n! / I > 
where Z denotes the set of integers. It is tempting to conjecture that the sets with basic represen- 
tations in an effective topology {B, v) are exactly the Gg sets, but there is a trivial counterexample 
to this conjecture if the effective topology is not Ti. 
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Theorem 12.2. Let A he a set in an effective topology {B,i>). 

(1) If A is an open set, then A has a basic representation. 

(2) If A is a closed set, then A has a basic representation. 

Proof. Suppose that A is an open set in the effective topology {B, v) and let 

i? = { r e domg V \ v(r) C A } 

Clearly, R C donig v and ii x ^ A then there does not exist a local basis for x 
with Cx Q v{R)^ since no member of ^{R) contains x. Alternatively, \i x £ A then, 
by the definition of a basis, for each basis element Bi that contains x there exists 
some basis element B2 such that x d B2 Q Bi A. That is, li x E A then for each 
basis element Bi with x ^ Bi, there exists a basis element B2^{B£B\x£BC 
A } with B2 C Bi. It immediately follows that Cx = { B G B \ x £ B C A} is a 
local basis for x and Cx ^ i^{R). By definition, i? is a basic representation of A. 
Now, if A is a closed set in {B, v), then let 

R = {r £ dome v\ Af] v{r) 7^ } 

Clearly, R C domg v and li x £ A then there exists a local basis Cx for a; with 
Cx C Namely, Cx is the set of all basis elements that contain x. Alternatively, 

if 2; ^ A, then since A is closed, every local basis Cx for x contains a basis element 
that does not intersect A. Therefore, Cx ^ v{R). We may conclude, by definition, 
that i? is a basic representation of A. □ 

Although we use a basic representation R to represent a set of points in an 
effective topology, the following theorem demonstrates that there is, in general, no 
effective procedure (relative to R) for finding oracles for those points. Nevertheless, 
if we restrict our attention to certain special classes of basic representations i?, then 
effective procedures do exist. See Section [TSl 

Theorem 12.3. Let B he a countable basis for the standard topology ofK" and let 
V be a coding for the basis. Then there does not exist a partial recursive function 4> 
satisfying the condition that for every singleton set {x} C M" and for every basic 
representation Rx of{x} in the effective topology {B, v), the function Xm\(f){Rx, m)] 
is an oracle for x in {B, v). 

Proof. Let {B, v) be an effective topology as in the statement of the theorem and 
suppose, as an assumption to be shown contradictory, that there exists a partial 
recursive function </) satisfying the condition that for every singleton set {x} C R" 
and for every basic representation Rx of {x} in the effective topology (S, i/), the 
function \m\^(f>{Rx^ m)] is an oracle for x in (S, v). Now consider any two distinct 
points X £ M" and y £ M", and let Rx be a basic representation of {a;} in {B,v). 
Since the standard topology of R" is Ti , there must exist a non-negative integer k 
such that 4){Rx, k) is defined and 

y^iy{cj){Rx,k)) 

Next, choose a program for computing (/). Note that since the computation for 
4>{Rx, k) has only finitely many steps, only finitely many non- negative integers are 
tested for membership in Rx during the course of the computation. Let C be the 
collection of all i G N such that i £ Rx and such that i is tested for membership 
in Rx during the course of the computation of (j){Rx,k). Similarly, let D be the 
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collection of alH G N such that i ^ Rx and such that i is tested for membership in 
Rx during the course of the computation of (t>{Rx, k). 

Now, choose any oracle -0 for y in [B, v). Note that -(/"(N) is a basic representation 
for {y} in [B^ v). And because C and D are finite sets, 

Ry ^ (V'(N) U C) - D 

is also a basic representation for {y} in (B, v). It follows that (t){Rx, k) = (f){Ry, k), 
because whenever i is tested for membership in Rx during the course of the com- 
putation of 4i{Rx, k), i £ Rx if and only if i e Ry. Therefore, 

y^,y{^iRx,k))=iy{^{Ry,k)) 

But by the definition of (j), Am[0(i?j,, m)] is an oracle for y. Hence, 

y e iy{(t){Ry,k)) 

This is a contradiction, so the assumption must be false. The partial recursive 
function does not exist. □ 

13. Basic Representations of Physical Models 

Let R be the set of all real numbers. For any two sets A and B, let A x B = 
{ (a, 6) I a G A & 6 S i? } be the Cartesian product of A with B. We write A'' to 
denote the set formed by taking the Cartesian product of A with itself k many times. 
For example, A'^ = {A x A) x A. As with Cantor's pairing function, (a, b, c) is an 
abbreviation for ((a, b), c), and so on. Similarly, we define the Cartesian projection 
function w''^ so that vd'^[x\^ X2, . . . , = Xi for each positive integer i < n. 

A physical model {S, A) with finitely many observable quantities is said to be in 
normal form if and only if S" C M" and A = {w", tUj , . . . , '^n}- 

Theorem 13.1. The following two conditions hold for any physical model {S,A) 
with finitely many observable quantities. 

(1) {S,A) is observationally equivalent to a physical model in normal form. 

(2) (S, A) is isomorphic to a physical model in normal form if and only if {S, A) 
is a reduced physical model. 

Proof. Begin by noting that if (T, B) is a physical model in normal form, and if 
ti ^ t2 for any ti £ T and t2 £ T, then wY{ti) ^ ■ujf{t2) for some positive integer 
i < n. Therefore, by definition, every physical model in normal form is a reduced 
physical model. It immediately follows from Theorem II . 71 that if a physical model 
(5', A) is isomorphic to a physical model in normal form, then (S, A) is a reduced 
physical model. 

To prove condition [TJ suppose that {S, A) is a physical model such that A — 
{ai,a2, . . . , an}. Define 

T = { (ai(s),a2(s),...,a„(s)) | s G S' } 

and let B = {wi,vu2, . . ■ Note that {T,B) is a physical model in normal 

form. Also note that the function (p : S ^ T given by 

(/)(s) {ai[s),a2[s), . . . ,a„(s)) 

is a surjection, and that ai{s) = n7"(0(s)) for all s G S" and all positive integers 
i < n. Therefore, there is an epimorphism from {S, A) to the reduced physical 
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model (T, B). Trivially, there is also an epimorphism from (T, B) to itself. We may 
conclude that {S,A) is observationally equivalent to {T,B). 

To prove condition[51 consider the special case where (S", A) is a reduced physical 
model. Because [S, A) is reduced, we have that if si ^ S2 for any si € S* and S2 € S, 
then there exists a positive integer i < n such that ai{si) 7^ ai(s2)- This implies 
that if si 7^ S2 then 0(si) 7^ 0(s2)- Hence, is an injection. Since (p is also a 
surjection, is a bijection. Therefore, if (S, A) is a reduced physical model, then 
{S, A) and (T, i?) are isomorphic. We have already proved the converse, that if 
{S, A) is isomorphic to a physical model in normal form, then (S, A) is a reduced 
physical model. Hence, condition [2] holds. □ 

The notion of a basic representation of a set can be generalized so that we may 
speak of basic representations of physical models in normal form. Given a physical 
model {S, A) in normal form with A = {ro", W2, . . . , w"}, we may choose sets Xi, 
X2, . . . , Xn such that wf{S) C Xi CM. for each positive integer i < n, and we may 
impose effective topologies {Bi, i^i), {B2, ^'2), • ■ • , {Bn, i^n) on these sets0 Define 

{Bi,i^i)(E}{B2,iy2)<»---(E){B 

to be the effective topology with basis B such that 

S = { Bi X S2 X • • • X B„ I Bi e Si & B2 e ^2 & • • • & B„ G S„ } 

and with coding ly such that 

iy(ai,a2, ■ ■ ■ ,a„) = 1^1 (oi) x 1/2(02) x • • • x i/„(a„) 

We call {B, u) the effective product of {Bi, vi), {B2, 1^2), ■ ■ ■ , {Bn, Vn)- Note that S is 
a set of points in the effective topology (S, v). A physical model (i?, H) is said to be 
a basic representation of the physical model (S", A) if i? is a basic representation of 
S in the effective topology (S, v) and if H is the set {tt", ttj , . . . , tt^} of projection 
functions. 

For example, the non-discrete continuous computable physical model of plane- 
tary motion (Model 14. ip is a basic representation of the simple model of planetary 
motion f Model [TTT|) . In particular. Model HTTl is obtained by imposing the effective 
topology (Iio,c , i) on the time in Model [TTTl where c — and by imposing the 
effective topology described in Footnote S] on the angular position in Model 11.11 
The product of these topologies is the topology for the surface of a cylinder. The 
states of Model [TTT] are a spiral path on the surface of that cylinder, and the set of 
states of Model 14.11 is a basic representation of the path. 

14. Data and Predictions 

In order to make predictions, we are often interested in finding the set of all 
states of a physical model which could account for a given collection of simultaneous 
measurements. That is, given a physical model {S,A) with A = {ai,a2,ot^, ■ ■ .}, 
and given real numbers xi, X2, ... , Xk^ we are interested in the set 

P ~ { s <^ S \ ai{s) = xi <k 0:2(5) = X2 ■ ■ ■ & ak{s) = Xk} 

^^If {S, A) is faithful, then the bases for these topologies are uniquely determined by the phys- 
ical operations used to measure each of the observable quantities. For example, if an observable 
quantity is an angle measurement, then the corresponding topology is the topology of a circle, and 
each basis element corresponds to a particular reading on the instrument that is used to measure 
angles. The idea that basis elements correspond to the values of measurements appears to have 
originated with reference I27| . 
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In this context, the real numbers xi, X2, ■ ■ ■ , Xk arc said to be the data, and P is 
the corresponding set of states predicted by the model. 

The following theorem shows that if we are given a basic representation of a 
physical model (5*, A) in normal form, together with complete oracles for the real 
numbers xi, X2, ... , Xk, then there is an effective procedure for finding a basic 
representation of the set P, provided that the underlying topology is Ti. (This is a 
rather weak requirement, since almost all topologies with practical applications in 
the sciences are Ti.) 

Theorem 14.1. Let {S, A) be a physical model in normal form with A = {zui, , 
. . . , roJJ} and let {R, H) he a basic representation of {S, A) in a T\ effective topology 

{B, v) = (Bi, ZVi) ® (S2, 1/2) O • • • O {Bn, z^n) 

If k < n and ipi, (j)2, ■ ■ ■ , ipk are complete oracles for Xi, X2, ... , x^ in the effective 
topologies {Bi,vi), {B2,i'2), ■■■ , {Bk,Vk), then there is a basic representation of 

P = {seS\ zu'^is) = xik m^{s) ^X2k ■■■ & tz7^(s) = Xk } 

in {B,h') that is recursively enumerable relative to R, (pi, (f)2, ■■ ■ , (j^k uniformly. 

Proof. Let the variables be defined as in the statement of the theorem and note 
that the set 

Q = { r e i? I (Vi e {1, 2, . . . , fc}) (3 m € N) [n^r] = ^m)] } 

is recursively enumerable relative to i?, (pi, (f)2, ... , (pk uniformly. (In fact, Q 
is recursively enumerable relative to ^, (j)2, ... , (pk uniformly, where ^ is a 
function that merely enumerates the members of R.) We claim that Q is a basic 
representation of P in {B, v). Since Q ^RC. domg v, it suffices to prove that s G P 
if and only if there exists a local basis Cg for s with Cg C v{Q). Or equivalently, 
it suffices to prove that s S P if and only if there exists an oracle ip for s with 
V'(N) C Q. 

Suppose s G P. Because i? is a basic representation of S in {B,v), there is an 
oracle for s in {B, u) such that tfi{H) C R. Moreover, for each positive integer 
i < k, the set t'j(7rf(V'(N))) is a local basis for wf{s) = Xi. And since (pi is a 
complete oracle for Xj, we have that 7r7(i/'(N)) C (piiN). Hence, if r = ipil) for 
some I e N, then there exists an m € N such that 7r"(r) = (pi{m). Therefore, by 
the definition of Q, %p{N) C Q. 

Conversely, suppose that tp is an oracle for some point s in (B,!/), and that 
V'(N) C Q. Then, for each positive integer i < k, we have that 7r"(^/'(N)) C (pi{N). 
Of course, i'i (7r"(V'(N))) is a local basis for wf{s) because i'[ip{N)) is a local basis 
for s. And by the definition of (pi, j/i((/)i(N)) is a local basis for Xi. Hence, a local 
basis for ro"(s) is a subset of a local basis for Xi in the effective topology {Bi, Vi). 
But because {B, v) is a Ti effective topology, [Bi, Vi) is also Ti. In a Ti topology, 
local bases for any two distinct points zi and Z2 must contain basis elements B\ 
and B2, respectively, such that Z2 4- ^1 ^"^^ zi ^ B2. Therefore, since a local basis 
for zuf{s) is a subset of a local basis for Xi, it must be the case that a7"(s) = Xi. 
We may conclude, by the definition of P, that s G P. □ 

A set 5 C M" is said to be the graph of a function : M'^ if and only if 

S = {{xi,...,Xk,Xk+l,...,Xn) eM" I tp{xi,...,Xk) = {Xk+l,...,Xn)} 
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And we say that a physical model (5, A) is induced by a function ip : R'^ M"^'^ if 
and only if (S, A) is in normal form and S is the graph of "0. Therefore, if {S, A) is 
induced by V' : — > R""'^ and we are given k real numbers xi, X2, ■ ■ ■ , Xk as data, 
then the corresponding set of states predicted by {S, A) is a singleton set P = {s}. 
Namely, 

S = {xi,X2, ■ ■ ■ , Xk,Xk+l,- ■ ■ ,Xn) 

where Xk+i, Xk+2, ■ ■ ■ ,Xn are the real numbers uniquely determined by the equation 

ip{xi,X2,...,Xk) = (a;fc+i,Xfc+2, . . . ,x„) 

It then follows from Theorem 114.11 that given complete oracles for xi^ X2, ■ ■ ■ ^ X}~ 
in the standard topology of R, and given a basic representation R of (S', A) in the 
standard topology of M" , there is an effective procedure (relative to the given oracles 
and R) for finding a basic representation of {s}. But by Theorem 112.31 there is, in 
general, no effective procedure for finding an oracle for s. In the next section we 
describe a special class of basic representations for which such an effective procedure 
does exist. 

15. Kreisel's Criterion 

A common way to interpret Kreisel's criterion is to say that a physical model 
{S, A) satisfies Kreisel's criterion on R'^ if and only if (S*, A) is induced by a function 
■0 : R'^ R"^*"' for some positive integer n > k, and for each positive integer 
j < n — k there is a partial recursive function Kj such that if ^i, 4>2, ■ ■ ■ , (f'k are 
nested oracles for real numbers xi, X2, ■ ■ ■ , Xk iri the effective topology (Z, t), then 
Kj((/>i, 02, . . . , (/>fc, m) is defined for all m e N and Am[Kj(0i, 02, ■ • ■ , 0fc, ™)] is a 
nested oracle for zz7"~'''(-0(xi, 2:2, . . . , Xk)) in (I, t). Note by Theorem II 1 .41 that the 
effective topology (I, t) in this statement can be replaced, without loss of generality, 
with any effective topology {B, l) that has a recursively enumerable subset relation 
and such that S C I is a basis for the standard topology of the real numbers. 

Practical computer models that use multiple-precision interval arithmetic |15j 
provide examples of physical models satisfying Kreisel's criterion. Typically, such 
models are induced by a function ip : 'MJ' R""*^ where, for each positive integer 
j < n — k, there is a recursive function such that if the data xi, X2, ... , Xk lie 
within the intervals (ai ; &i), (02 ; 62), • ■ • , (ofe ; bk) respectively, then ■n?""'" ('0(a;i, X2, 
. . . ,Xk)) lies within the interval 

^j{{ai;bi),{a2;b2),.-.,{ak;bk)) 

In such a case, the partial recursive function Hj in Kreisel's criterion is given by 

Kj(0i,02, ■.■,(t)k,m) = ^j(0i(TO),02(m), . . .,(t)k{m)) 

We are now prepared to state the following theorem, which holds uniformly. 

Theorem 15.1. If a physical model satisfies Kreisel's criterion on R'^, then the 
model has a basic representation that is isomorphic to a computable physical model. 

Proof. Suppose that {S, A) is a physical model satisfying Kreisel's criterion on R*^. 
In particular, suppose that {S,A) is induced by a function ij; : R'' ^ R"^'^, and 
for each positive integer j < n — k suppose there is a partial recursive function Kj 
such that if 01, 02, • ■ • , 0fc are nested oracles for real numbers xi, X2, ... , Xk in 
the effective topology (Iio,c , t), then Kj(0i, 02, . . . , 0fc, m) is defined for all m G N 
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and Am[Kj (01, 02, ■ • • , ™)] is a nested oracle for w"~'' (^'ip{xi,X2, ■ ■ ■ , Xk)) in 
(Xio,c , where c is a positive rational number. 

Let / = domxio c Then, for each interval u € I with midpoint p and with d 
digits of precision, define the partial recursive function ct„ so that 




Op{l) if / < d 

undefined if / > d 



for each / G N, where Op is the standard decimal oracle described in SectionH) Note 
that for each a; e M and each m e N, CTq^ (,„)(/) — Ox{l) for all non-negative integers 
I < m. Now, for each positive integer j < n — k, choose a program to compute 
Kj and let Kj (ctuj , cr„2 , ■ • ■ , cr^^. , m) be undefined if for some positive integer i < k 
and some I e N the program calls cr^.[l) in the course of the computation and 
CFui(l) is undefined. Note that given ui, U2, ■ ■ ■ , Wfc, the set of all m € N such 
that «;j((Tuj,(T„2i • ■ • i^'iikJ™) is defined is a recursively enumerable set, since for 
each m we can follow the computation and test whether or not (7^.(1) is defined 
whenever a;^.{l) is called by the program, for any i and /. Let R be the set of 
all U2, ■ ■ ■ , Un) such that Ui € I for each positive integer i < k, and such that 
Uk+j = Kj {<7ui , cTtia I • ■ • I o'tifc I ™) for some m e N if j < n — fc is a positive integer. 
Note that R is also recursively enumerable. Let H = {tt", , . . . , tt"}. We claim 
that {R,H) is a basic representation of {S,A). 

As a brief digression from the proof, suppose that (ui, U2, . . . , m„) G R and note 
that for each positive integer i < fc, if a;^ S i{ui) then there exists a nested oracle (f>i 
for Xi in (Xio,c , t) such that 4'i{m) — cr„. (m) for all m S N where CTu. (m) is defined. 
And since for each positive integer j < n—k we have that Xm [kj ((/>i ,(f)2,-.-,(f>k, m)] 
is a nested oracle for Xk+j — ^2, ■ ■ ■ , ^k)) , it follows that 

Xk+j e t(Kj(0i, 02, • • • , (j^k^m)) = t(Kj(cr„i,CT„2, . . . , au^,m)) = /-(w^+j) 

for some to G N. Hence, given any (wi, M2, . . . , Wn) G i?, if e for each 

positive integer i < k, then Xfe+j = (jp{xi,X2, ■ ■ ■ ,Xk)) G L{uk+j) for each 
positive integer j < n — k. 

Now, returning to the proof of Theorem 115.11 let 

n factors 

and note that the basis of (S, v) is a basis for the standard topology of M". Note 
that because (S', A) satisfies Kreisel's criterion, the function -0 which induces (5, A) 
is continuous. And since continuous real functions have closed graphs, the set S 
is closed in the standard topology of M". Because 5 is closed, to prove that R is 
a basic representation of 5 in (S, v)^ it suffices to show that for each r £ R there 
exists an a; G 5 with x G v{r), and that for each x £ S there is a local basis Cx for 
X with Cx Q v(R). 

By the definition of i?, if r = (iti, M2, . . . , u„) G R then for each positive integer 

j < n - fc, 

Uk+j = Kj{aui , cr„2, . . . , (T„j^, m) 

for some to G N. So, if pi is the midpoint of the interval Ui for each positive integer 
i < k , then 

Wfc+j = Kj (opj , , . . . , Opj, , to) 
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And by Kreisel's criterion Am [kj (op^ , Op^ , • • • , Op^ , m)] is an oracle for zu" ^ {ip{pi , 
P2, ■ ■ ■ ,Pk)) ■ Therefore, for each r G R there exists an {xi,X2, ■ ■ ■ ,Xn) S S with 
(xi,X2, ■ ■ ■ ,Xn) € t^l?")- Namely, Xi — pi for each positive integer i < k and 

Xk+j = zu'J^'' {t(;{pi,p2, . . . ,Pk)) 

for each positive integer j < n — k. 

Now suppose that (xi,X2, ■ ■ ■ ,Xn) is an arbitrary member of S. Again, by 
Kreisel's criterion, for each positive integer j < n — k, the function Am [kj {ox^ , , 
. . . , Oxf. , m)] is an oracle for Xk+j ■ But for each m e N and each positive integer 
i < k, the computation for Kj{oxi,Ox2, ■ ■ ■ TOxi^,m) has only finitely many steps, 
and so the oracle Ox^ can only be called finitely many times during the course of 
the computation. Hence, for each m e N there exists a non-negative integer li 
for each i < k, such that for any non-negative integer l[ > li, if Ui = 0x^(1'^) then 
Kj (au-i , cr„2 , . . . , cr„fc , m) is defined and 

Of course, for each positive integer j < n — k the interval 

Uk+J = (Oxi ,0x2: ■ ■ ■ :Oxk,m) 

can be made arbitrarily small by choosing a suitably large value of m, and for each 
positive integer i < k the interval Ui can be made arbitrarily small by choosing 
a suitably large value of l[. Furthermore, by definition, it2, . . . , m„) S R. It 
immediately follows that for each x = {xi,X2, . . . ,Xn) S •S' there is a local basis Cx 
for X such that Cx ^ i^(-R). We may conclude that {R, H) is a basic representation 
of (S, A). And since i? is recursively enumerable, it follows from Theorem 1 1 . 5 1 that 
(i?, H) is isomorphic to a computable physical model. □ 

A physical model [S, A) that satisfies Kreisel's criterion on R'^ is uniquely de- 
termined by the functions ki, k.2, ... , Kn-k- Moreover, the proof of Theorem 1 15. II 
describes an effective procedure for finding a basic representation of {S,A), given 
programs for computing ni, K2, ... , Kn-k- Let lCk,n,c be the collection of all basic 
representations of physical models that are constructed from physical models sat- 
isfying Kreisel's criterion according to the procedure in the proof of Theorem 115.11 
where c is the positive rational number which appears in that proof. An immediate 
question is whether there exists an effective procedure for the inverse operation. 
That is, given a basic representation in ICk,n,c, is there an effective procedure for 
constructing partial recursive functions ki, K2, ... , Kn-k'^ In the proof of the fol- 
lowing theorem, we show that the answer is "Yes." Therefore, for every physical 
model satisfying Kreisel's criterion on R*^, there is a computable physical model 
that may be used in its place, to predict the values of observable quantities given 
the data. 

Theorem 15.2. // (R, H) e lCk,n,c and if tp ■ ^ R"^*^ is the function whose 
graph has basic representation R, then there exist partial recursive functions ki, 

K2, ... , Kn—k 

such that if (j)2, ■ ■ ■ , (^k are nested oracles for real numbers x\ , 
X2, ■ . ■ , Xk in the effective topology (Iio,c , '•), then for each positive integer j < n—k, 
Kj{(f>i, (f>2, ■ ■ ■ ,(f>k,m-) is defined for all m € N and Am[Kj(0i, 02, • ■ • , </'/c, w)] *s a 
nested oracle for tu"~'^(V'(xi, 0:2, . . . , Xk)) in (Iio.c , t)- 
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Proof. Suppose that {R, H) E ICk.n.c and that i/; : R'^ — > R"^''' is the function whose 
graph has basic representation R. Note by the proof of Theorem 115.11 that R is 
recursively enumerable. Now, given any oracles 0i, 02, • • ■ , 0fc for real numbers 
xi, X2, ■ ■ ■ , Xk (Xio.c , i), it follows from Theorem 111.21 that there are complete 
oracles 4>[, 4''2, ■ ■ ■ , 4>'k for xi, X2, ■ ■ ■ , Xk va (Iio.c , '-), such that 0']^ , 02 1 ■ ■ ■ , 4>'k 
recursive relative to 0i, 02, ■ ■ • , ^fe uniformly. Let Xfe+i, Xk+2^ ■ ■ ■ , a^n be the real 
numbers uniquely determined by the equation 

lj}{xi,X2, ... ,Xk) = {xk+l,Xk+2, ■ . ■ ,Xn) 

Then by the proof of Theorem 114.11 

Q = { (ui , U2, . . . , Un) e i? I (Vi e {1, 2, . . . , fc}) (3 m e N) [u, = 0^ (m)] } 
is a basic representation of {{xi,X2, . . . ,Xn)} in 

n factors 

{B, v) = (Xio^c , ® (Iio.c i^W^c , 

And since R is recursively enumerable, the set Q is recursively enumerable relative 

01, 02, . . . , 0fc uniformly. 

Now, since Xi e L{(f)i(rn)) for each positive integer z < fc, it follows from the 
definition of Q that Xi € for each (mi,U2, . . . ,Un) € Q- But recall from the 

proof of Theorem 1 1 5 . 1 1 that R has the property that if xi e for each positive 
integer i < k, then G iiuk+j) for each positive integer j < n — k. Hence, 

(xi,a;2, . . . ,x„) e iy{ui,U2,...,Un) 

for each (ui, U2, ■ • ■ , w„) G Q. It immediately follows from the definition of a basic 
representation that v{Q) is a local basis for the point (a;i,a;2, . . . ,Xn). Therefore, 
for any function k : N — i? such that k(N) — Q, the function k is an oracle for 
{xi,X2, . . . ,Xn) in {B,!^). And since Q is recursively enumerable relative to 0i, 

02 , ... , 0fe uniformly, there is a nested oracle k that is recursive relative to 0i , 
02, . . . , 0fe uniformly. So, if we define 

Kj(0i,02, . . . ,0fe,m) = 7r^+j(K(m)) 

for each positive integer j < n — k and for each m G N, then Kj is partial recursive, 
Kj (01 , 02 , . . . , 0fc , m) is defined for all m G N, and Am [nj (0i , 02 , . . . , 0fe , to)] is a 
nested oracle for n7"~'^(-0(a;i,a;2, ... jXfc)) — Xk+j in (Iio,c , '-)• D 
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